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Abstract 

Gamma-ray bursts are generally considered to be the result of internal shocks 
generated in an inhomogeneous relativistic outflow that arises from a fireball. In 
such shocks, the Fermi acceleration of protons is naturally expected to be at least as 
efficient as that of electrons. We investigate the consequences of proton acceleration 
in the residual thermal photon field of a fireball, especially those on the emission 
spectra of photons. In contrast to most other studies, we do not invoke a direct 
electron acceleration in shock waves. We show that the residual radiation field of the 
fireball can ignite the photopion production and subsequent cascade, and that the 
photons emitted in this process further enhance the photon-initiated cascade. We 
find that Fermi accelerated protons with > 10 13 eV efficiently bring their energy into 
pion production and subsequent photon emission. The particles cascading from the 
pions emit photons over a wide energy range. The photons emitted from electron- 
positron pairs distribute continuously from the GeV range down to the X-ray range, 
while muons can emit gamma-rays by synchrotron radiation with a break at around 
1-10 MeV. We also discuss several radiation processes which may possibly produce 
a break feature in the MeV range, as observed in addition to muon synchrotron 
radiation. 

Key words: gamma rays: bursts — radiation mechanisms: non-thermal — shock 
waves 



1. Introduction 



A widely accepted scenario for producing gamma-ray bursts (GRBs) is dissipation of the 
kinetic energy of a relativistic flow by relativistic shocks (see, e.g., a review by Piran 1999). The 
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rapid time variabilities observed require that the GRB itself must arise from internal shocks 
within the flow, while the afterglow is due to external shock produced as the flow is decelerated 
upon collisions with the ambient medium. Although how the central engine of GRB produces 
the relativistic flow remains an open problem, the standard scenario is the high-entropy fireball 
scenario (Shemi, Piran 1990; Rees, Meszaros 1992; Meszaros, Rees 1993), in which a large 
amount of internal energy is injected into a compact region, such that the internal energy per 
one proton far exceeds the rest mass of the proton (Goodman 1986; Paczyhski 1986). The huge 
internal energy of the fireball will be converted into the bulk kinetic energy of the relativistic 
flow, which we call the fireball wind. This fireball wind is accelerated to a relativistic velocity 
and the internal energy of the wind is finally dominated by baryonic matter. It is assumed that 
the fireball wind is inhomogeneous in order to produce internal shocks. The inhomogeneity in 
the wind may grow, and it ends up with multiple shells with different bulk Lorentz factors. 

It is to be noted that protons carry a much larger amount of energy than electrons 
in the flow. In the standard GRB model, the observed spectra of GRBs are interpreted in 
terms of the synchrotron radiation emitted by relativistic electrons, which requires that a large 
fraction of the kinetic energy carried by protons is efficiently converted into that of relativistic 
electrons in the shocked region. However, this premise has not been theoretically proven. It is 
apparent that the Coulomb interaction cannot transport the internal energy of heated protons 
into electrons to achieve energy equipartition (Totani 1999), because the time scale of the 
Coulomb interaction is much longer than the dynamical time scale. As long as protons and 
electrons are accelerated independently in the shocked region, the total energy of accelerated 
protons should overweigh that of accelerated electrons by an order of their mass ratio. Although 
more efficient electron acceleration mechanisms in shocks may exist, at present no convincing 
one has been specified. On the other hand, efficient proton acceleration may lead to high- 
energy photon emission through photon-initiated cascades, as investigated in this paper. Even 
if we regard that efficient electron acceleration by shocks is quite possible, proton acceleration 
is naturally expected to be at least as efficient, and most probably more efficient; thus, the 
consequences of proton acceleration on photon spectra should be investigated. 

One of the most important characteristics of GRB spectra is the existence of a typical 
break energy scale. The observed spectra of GRBs are approximated by a broken power law 
and the photon number spectra are approximated as ~ e~ 2 above the break energy and ~ e~ l 
below that. The apparent clustering of the break energy of GRB spectra in the 50 keV-1 MeV 
range (Preece et al. 2000) is an important key to understanding the physical processes in GRBs. 
Within the scheme of the standard model, the clustering of the break energy and the paucity 
of low energy photons below the break in GRB spectra (Preece et al. 1998, 2000) remain to be 
explained. The break energy scale has been conventionally adjusted by choosing the minimum 
energy of accelerated electrons. It is thus interesting to see if proton acceleration can provide 
other mechanisms to reproduce the break energy scale. 
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Although direct synchrotron emission from protons is too weak and too soft to reproduce 
GRBs for the expected values of the magnetic field, the emission from secondary particles is not 
safely neglected, provided that the Fermi-accelerated protons account for a large fraction of the 
internal energy and that there exist copious target photons to trigger the cascades. In order to 
optimistically estimate the high-energy neutrinos and extremely high-energy cosmic rays from 
GRBs, the Fermi-accelerated protons have often been assumed to have a large fraction of the 
internal energy. The physical conditions in GRBs imply that protons may be Fermi-accelerated 
to high energies (Waxman 1995). Waxman, Bahcall (1997) pointed out that a large fraction 
of the fireball energy is expected to be converted into high-energy neutrinos by photomeson 
production from the accelerated particles. In this paper, we explore the consequences of proton 
acceleration on the photon spectra and investigate if GRBs themselves are produced by proton 
acceleration. 

Assuming that the accelerated protons account for a large fraction of the internal energy, 
we discuss the emission from particles produced in the cascades from high-energy protons. For 
simplicity and clarity, we assume that the energy density of directly accelerated electrons is 
very small compared to that of protons, so that the emission from this component is negligible 
in contrast to that found in almost all previous studies. In section 2, we argue that the 
residual radiation field of the fireball can trigger photopion production of relativistic protons. In 
section 3, the emission from the cascading particles is investigated. In section 4, we discuss the 
possibility to explain the observed spectral property of GRBs. Our conclusions are summarized 
in section 5. 

2. Cascade Initiated by a Residual Radiation Field of the Fireball 

In this section we discuss photopion production in the internal shock against residual 
thermal photons of the fireball. While most of previous studies have ignored the effects of 
the residual photons, we argue that this is not legitimate, and that these photons play an 
important role to trigger photopion production. The cascade triggered by the collisions of 
shock-accelerated protons with residual photons can efficiently transport the dissipated energy 
into electron-positron pairs. A schematic diagram of the energy transfer is shown in figure 1. 
In the following, we discuss the ignition of GRBs according to this diagram. 

In order to discuss this issue, we need to specify the distribution of the residual radiation 
field. We regard that the central engine produces an inhomogeneous fireball wind for a long 
duration compared to the dynamical time scale of the central engine. The inner region of 
the fireball wind is optically thick and the radiation and matter in the fireball wind behaves 
like a single fluid moving with the same velocity. In the radiation-dominated stages, as the 
wind expands, the Lorentz factor of the wind increases in proportion to the radius (Piran et 
al. 1993). After the fireball wind becomes matter-dominated, the wind plasma coasts with 
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Fig. 1. Process of the cascade ignited by the residual radiation in the fireball (see section 2). We neglect 
the effects of direct electron acceleration. The solid arrows show the energy flow. The dashed arrows 
indicate the interactions that induce energy transfer to lower energy particles. The typical energies of each 
particle measured in the comoving frame are noted. 

a constant radial speed. The fireball wind becomes optically thin at the photospheric radius 
where matter and radiation decouple. We treat this decoupled radiation field as the residual 
radiation field. 

In the internal shock scenario, there exist inhomogeneities of the matter in the wind, 
and they end up with multiple shells with different bulk Lorentz factors. Internal shocks arise 
from collisions of the shells. As for the Lorentz factor of the shells, we follow the general 
constraints obtained in the conventional internal-shock scenario. The bulk Lorentz factor of 
the shells lies between one hundred and a few thousand (Meszaros, Rees 2000). Moreover, the 
energy-efficiency argument concerning the internal shock model requires a large dispersion of 
the Lorentz factor of individual shells (Kobayashi, Sari 2001; Guetta et al. 2001). Efficient 
dissipation of the kinetic energy is possible only when the difference in the Lorentz factors of 
the two colliding shells is sufficiently large. Thus, we regard that independent multiple shells 
are arising intermittently from the central engine to satisfy the required large dispersion of their 
Lorentz factors. 

On the other hand, the inhomogeneity of the residual radiation field is less clear and 
the degree of inhomogeneities largely depends on the character of the central engine, about 
which we have little information. The residual photon field in the wind may not necessarily 
distribute discontinuously as do the matter shells. Therefore, we consider two extreme cases for 
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Fig. 2. Schematic representation for the Discrete Model (above) and the Continuous Model (below). The 
white shell (above) and the dashed lines (below) represent fluid outside of i? p h- 

the residual photon field (see figure 2). In one case, we assume that radiation is strongly coupled 
with matter and is contained in the shell. In this case, each shell is regarded to represent an 
independent fireball. We assume that there is no matter and radiation between the shells before 
the radiation field in a shell decouples with matter. In this case, we completely neglect the 
radiation field emitted from other shells. Hereafter, we call this case the "Discrete Model" . 

In the other extreme case, we assume that the photon field is given by a quasi-steady 
continuous wind, for which the effects of the inhomogeneities are neglected. Inside the pho- 
tosphere the radiation and matter are considered to be continuously distributed and coupled 
with each other. Outside of the photosphere photons escape freely to infinity. In this case the 
radiation field behaves smoothly, as determined from the wind parameters. Even though this 
treatment formally neglects the inhomogeneities in the wind, in reality inhomogeneities should 
exist, and we assume that the internal shocks occur and high-energy protons are injected at 
some places in the radiation field. We do not discuss the detail of the internal shocks in this 
model. The latter case is called the "Continuous Model" hereafter. 

For simplicity, we assume spherical symmetry for both models. 

2.1. Residual Photon Field 

2.1.1. Discrete Model 

We first discuss the Discrete Model. We consider the situation where the inner faster 
shell overtakes the outer slower one and generates forward and reverse shocks in the respective 
shells. We suppose that the typical Lorentz factors of the faster and slower shells are T r ~ 
1000 and T s ~ 100, respectively. Each shell is emitted from the central engine as a fireball with 
some initial radius A. Each fireball ends up as a coasting matter shell, retaining its width (A) 
in the observer frame. The radiation field in each shell evolves independently of other shells. 
In the initial radiation-dominated phase, the radiation energy density in the fireball behaves 
as e oc r -4 (r is the location of the shell), and the bulk Lorentz factor of the shell increases in 
proportion to r (Piran et al. 1993). After the shell enters the matter-dominated regime, as long 
as the shell is optically thick to scattering, the radiation energy density declines as e oc r -8 / 3 
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and the shell eventually coasts with T = E/Mc 2 , where M and E are the total baryon rest mass 
and the total energy of the shell, respectively. Since there is no matter between shells in the 
Discrete Model, once a photon escapes from the shell, it moves freely. 
In this case, the shell becomes optically thin at the radius 



where <tt is the Thomson cross section. The radius R p h becomes 6 x 10 12 -6 x 10 13 cm for a 
range Mc 2 = 10 48 -10 50 erg. The kinetic energy of one shell which we consider corresponds to 
10 50 -10 53 erg for T = 100-1000. Then, we obtain the temperature of the residual radiation field 
at r = R p h in the shell rest frame as 



where T 3 = T/10 3 , A 7 = A/10 7 cm, and M 48 = Mc 2 /10 48 erg. The parameter dependence of 
the temperature is very weak. The shell coasts R esc = 2T 2 A ~2x 10 13 A 7 r| cm from the radius 
r = i? p h before the residual photons escape from the shell. Thus, the energy density of residual 
photons in the comoving frame is as large as 10 10 erg cm -3 when the shell reaches a few times 
10 13 cm for T = 10 3 . 

The inner faster shell overtakes the outer slower one at r ~ T 2 5, where 5 is the initial 
separation between the two shells in the observer rest frame. When 5 is small enough, the 
collisions occur well inside of _R p h and no major observational feature is expected. For 5 ~ 10 9 
cm (which is consistent with the observed time variability of GRBs) and T s > 100, the two 
shells collide outside of i? p h ~ 10 13 cm. The typical radius where the internal shocks occur can 
be comparable to R p h, unless the initial separation of two shells is limited to larger than 10 10 
cm. Since the collision radius is comparable to R esc for the faster shell with T r ~ 10 3 , internal 
shocks can occur in the thermal radiation field in the faster shell with temperature T ~ 100 
eV. On the other hand, when 5 is much larger than 10 10 cm, internal shocks occur far outside 
of i? p h and the residual photons have already escaped from the shell. In this case there is no 
photon field in the shell when the internal shocks occur, which is not interesting for the present 
study. In the Discrete Model, therefore, we consider only the internal shocks that occur around 
the radius r = R p ^. In the comoving frame the width of the shocked region in the outer shell is 
thinner than that of the inner shell when A and M of the two shells are similar to each other. 
Since the photons in the shocked region in the outer shell escape more promptly, we discuss 
only the reverse shocked region of the inner shell. Hereafter, we denote r 3 as the normalized 
Lorentz factor of the inner shell (r r = 10 3 r 3 ). 

The internal shock is mildly relativistic, i.e., the Lorentz factor Tr of the reverse shock 
observed from the rest frame of the inner shell is a few. For example, when two shells with 
the same M and A collide with T s = 100 and T r = 1000, we have = 2.57. From the jump 
condition (Sari, Piran 1995), the shell width in the comoving frame of the shocked matter is 




(1) 



T = 130(r 3 A 7 M 48 ) 



1/12 



eV, 



(2) 
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/ = Ar r /(4rR + 3) — 10 9 Ayr3 cm. Then, the internal energy density in the shocked shell at 
r = i? p h is obtained as 

where we have assumed (Tr — 1)(4Fr + 3) ~ 10. Thus, the energy density of residual photons 
is about 1% of the internal energy of baryons. 

In the internal or external shocks of GRBs strong magnetic fields should be produced, 
since the expected magnetic field frozen from the central engine is several orders of magni- 
tude smaller than the value required to produce gamma-rays because of the adiabatic expan- 
sion. Although no specific mechanism is widely accepted, the relativistic two-stream instability 
(Medvedev, Loeb 1999) is one of the candidates to generate such a high magnetic field. In 
this paper, we assume that the constant ratio (/b) of the internal energy ({Tin) goes into the 
magnetic field; then, the magnetic field at radius i? p h becomes 

B(R ph ) = vWet/in ^ 2 x 10 6 (J^J ' Af 1/2 r 3 - 1/2 G. (4) 

2.1.2. Continuous Model 

Differently from the Discrete Model, for the Continuous Model we consider a quasi- 
steady fireball wind of total luminosity L = 10 52 L 52 erg s _1 expanding from the initial radius A. 
Although the radiation and matter are assumed to be continuously distributed in the estimation 
of the distribution of the physical quantities, in reality build-in inhomogeneities can make the 
internal shocks occur somewhere and high-energy protons are injected. We do not deal with the 
details of the internal shocks in this model. The injected high-energy protons will interact with 
the continuous radiation field. The wind is characterized by the average baryon load parameter, 
rj = L/Mc 2 = 100?72. The photosphere and photospheric temperatures in this case were well 
studied by Meszaros, Rees (2000). In the radiation-dominated region the bulk Lorentz factor 
of the wind increases as T oc r (r is the radial coordinate in this model). For relatively low 
values of n < 77* = (LaT/47rm p c 3 A) 1//4 = 10 3 L52 4 A 7 1//4 , the increase of T stops at the saturation 
radius, r s = 77 A, and the flow continues to coast with r = 77 above r s . In this paper we consider 
only this range of 77. The photosphere is located at a radius 

We obtain the radiation temperature at r = i? p h in the wind rest frame as 

T = 40t7 2 5/3 L 52 5/12 A^ /6 eV. (6) 

Although the representative numerical values of i? p h and T are similar to those of the Discrete 
Model, their parameter dependence is rather strong compared with the dependence in the 
Discrete Model. Outside of i? p h, the radiation field is approximated by a diluted black body 
with temperature T. In the wind rest frame, the number density of the residual photons declines 
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as r~ 2 . Although the energy density of residual photons is an order of magnitude smaller than 
that in the Discrete Model, the photon distribution is not limited in narrow shells, and they 
are available for a longer time scale. 

The average total energy density in the wind rest frame, U = Lj ' ^r 2 ifc^ is given by 



These values are for a continuous wind, and the effects of inhomogeneities are not taken into 
account. When an inhomogeneity causes an internal shock which propagates in the background 
wind flow, the internal energy and magnetic field can be larger than the values given in equations 
(8) and (9) locally. However, for simplicity, in the following numerical estimation we do not 
consider such factors in the Continuous Model. The dynamical time scale in the Continuous 
Model is taken to be r/cq, which is typically two orders of magnitude longer than that in the 
Discrete Model l/c. 

2.2. Photopion Production by Accelerated Protons 

The maximum energy of accelerated protons is estimated from the conditions that the 
cooling time of protons should be longer than the Fermi acceleration time scale, and that the 
Larmor radius should be smaller than the shell width. For the Discrete Model, the maximum 
energy of protons at r = i? p h is most likely determined by the latter condition when /b<1. 
For /b = 0.1 and I = 10 9 cm, however, both of the conditions give comparable values of the 
maximum energy, which turn out to be ~ 10 17 eV. For the Continuous Model, a similar value of 
the maximum energy is obtained. The accelerated protons would obey a power-law distribution 
and may account for a large fraction of the internal energy. The thermal photons interact with 
these protons to produce photopions. 

Although in the shocked fluid rest frame the residual photons undergo a red-/blue-shift 
and their distribution becomes anisotropic, an isotropic distribution is a quite good approxi- 
mation in a relativistic flow. Thus, we neglect the anisotropy and suppose that the accelerated 
protons interact with the isotropic thermal radiation field with a given temperature. The typ- 
ical Lorentz factor of accelerated protons that interact resonantly with the thermal photons 
is 7p,t yP = £th/T, where e t h — 145 MeV is the threshold energy of a photon in the proton rest 
frame. 




(7) 



~ 2 x 10 10 rilL^n- 2 erg cm" 3 , 
where 71 = r/_R ph . If the fraction f B of U goes into the magnetic field, it becomes 

ft \ V2 




(8) 



(9) 
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2.2.1. Discrete Model 

For the Discrete Model, we consider the internal shocks that occur around only the 
radius r = i? p h. The typical energy of pion producing protons becomes -Ep )typ — 10 15 eV for 
T = 100 eV. At r = i? p h the number density of the photons is written as 

S7T £ 2 

n( ^ )= (^exp(. 7 /T)-r (10) 
The time scale of photopion production, t n , is written as 

rn rco Tl(F \ 

t~\ lp ) = 2vrc / d0sin0(l-cos0) / de 1 --^a 7T (x), (11) 

J0 Je' th 47T 

where 9 is the incident angle between the two interacting particles. We approximate the 
photopion production cross section (e.g. Stecker 1968) by a broken power-law profile as 
a^ x ) = 5 x 10- 28 (x/590) 3 ' 2 cm 2 for 290 < x < 590 and cr n (x) = 5 x 10- 28 (x/590)-°' 7 cm 2 for 
590 < x < 9800, where x m cC 2 is the photon energy in the proton rest frame. In the fluid rest 
frame the energy of a photon is expressed as e 7 = x m cC 2 / 7 P (1 — cos0) . Since the residual photons 
escape in the dynamical time scale l/c, we regard that efficient pion production occurs if t n is 
less than l/c. In figure 3, we plot the time scale t n and = t~ x l/c for the photopion production 
in the 100 eV radiation field. It is seen that for E p = E ptyp ~ 10 15 eV T n becomes ~ 8(//10 9 cm), 
which means that efficient pion production occurs. The threshold energy of protons for pion 
production, E PtT , at which r n becomes one, is 2 x 10 14 eV for T = 100 eV. Therefore, protons 
with 10 14 -10 17 eV efficiently produce 7r°s and 7r ± s. The photopion production would not occur 
outside of a few times _R p h because the photon number density decreases as r n oc r~ 2 , even if 
the residual photon field remains in the shell. Since the average inelasticity is about 20%, the 
typical energy of produced pions should be 10 13 -10 16 eV. 

2.2.2. Continuous Model 

When an internal shock occurs in a continuous wind, the accelerated protons will be 
injected into the wind. The accelerated protons can interact with photons during much longer 
period of time (r/crj) than in the Discrete Model {l/c). In figure 4 we plot the threshold energy, 
E p ^ T , at the photosphere calculated numerically in a similar way to that in the Discrete Model. 
The energy E p T weakly depends on L 52 . The dotted line in figure 4 corresponds to 77 = 77*, 
the maximum value allowed for the saturated fireball wind. For fiducial values of 7/2 = 1 and 
L52 = 1, the threshold energy becomes 3 x 10 14 eV, which is almost the same as the typical value 
for the Discrete Model. Thus, in the Continuous Model, too, efficient photopion production 
occurs around the photosphere. 

This may be understood qualitatively as follows. Because r n is a steeply increasing 
function of the energy in the low-energy range, as can be seen in figure 3, the threshold energy 
is roughly taken as a fixed fraction of the typical energy for photopion production, which 
is inversely proportional to the photon temperature. Since r n (E Pjtyp ) oc T 3 ■ Rp h /r), we may 
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Fig. 3. Time scale (solid line) and r w (dotted line) of protons for photopion production in the thermal 
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Fig. 5. Threshold energy, E V>T , for r/2 = 1 and A7 = 1 against r in the Continuous Case. The dotted line 
is the same as in figure 4. 

approximate E PyT oc E PttJP oc L?^ 2 ■ It is seen that this simple relation well reproduces the 
plots in figure 4. 

Since in the Continuous Model the residual photons distribute continuously, the photo- 
pion production can occur even outside of i? p h. In figure 5, we show the threshold energy, E PjT , 
against r for 7/2 = 1. Although the value t w drops as oc r _1 , E PtT itself does not increase drasti- 
cally. The region where efficient photopion production can occur by the thermal radiation field 
is inside of r ~ 100i? p h. Therefore, we expect a large fraction of the energy of protons 10 14 
eV for L 52 ~ 1) to efficiently go to the energy of pions. 

In the Discrete Model we have assumed that the internal shocks occur near r = i? p h. 
Therefore, the physical conditions in the shocked region are somewhat restricted. On the other 
hand, in the Continuous Model, photopion production occurs over a wide range of r and the 
properties of the photon emission may depend on the parameters 77, L, and 1Z. For simplicity, 
we describe the emission from the cascading particles mainly along the Discrete Model below, 
while the Continuous Model is briefly mentioned. As will be shown later, however, there are 
no drastic differences in the results between the two Models. 

2. 3. Behavior of Pions 

The created 7r°s decay immediately into two gamma-rays, which are converted into 
electron-positron pairs, while the minor portion directly escapes from the shocked region. The 
fate of these high-energy photons is discussed in the next subsection. On the other hand, the 
lifetime of most 7r ± s (E n > 0.2E ptyp ~2x 10 14 eV) is longer than the dynamical time scale 
l/c ~ 0.03 s for the Discrete Model because of the large Lorentz factor. The energy density 
of residual photons is about 10 10 erg cm -3 , which is smaller than the energy density of the 
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magnetic field if B is larger than 5 x 10 5 G. Thus, synchrotron loss is the dominant cooling 
process for the charged pions. 

The synchrotron cooling time scale of 7r ± is expressed as 

Qirmlc 3 / / B V / E w \ ~ 



U,cooi = ~2xl0- d - " A 7 r 3 sec, (12) 

o r T m^i? 2 7 7r VO- 1 / V2xl0 14 eV/ 

where 7^ is the Lorentz factor of pions. The cooling time of high-energy pions becomes shorter 

than their lifetime. The lifetime of pions becomes comparable to the cooling time for 

2 _ formfe _ 2 

where t n>0 = 2.6 x 10 -8 s is the lifetime of pions in the rest frame. Thus, charged pions (E n ^ 
2 x 10 14 eV) cool down to m 7r c 2 7 7r o = 5 x 10 13 (/B/0.1) _1 ^ 2 eV before they decay into muons. 
Charged pions with lower energies decay into charged muons before significant cooling, and 
most of their energy is converted into photons by synchrotron radiation because the lifetime of 
charged muons is much longer than that of charged pions. 

The energy of the synchrotron photons emitted by charged pions is 

e =^ 7 2 ~9xW E " ) 2 (—) eV (14) 
* v2xl0 14 eW Uo 6 G7 1 ' 

Thus, pions of 10 14 -10 15 eV emit photons of 10 7 -10 9 eV. The energy of muon synchrotron 
radiation is distributed mainly in the range of 10 5 -10 7 eV. The conversion efficiency into high- 
energy photons in the cascade process is fairly large, much larger than that into neutrinos. 

In the Continuous Model, since the magnetic field is weaker the synchrotron cooling time 
scale becomes longer, comparable to the decay time for charged pions with energy of 10 15 eV. In 
a similar way to that in the Discrete Model we obtain m 7r c 2 7 7r0 — 5 x 10 14 (/b/0.1) _1 / 2 7/ 2 2 -^52 
eV, which is slightly larger than E P , T . However, the synchrotron cooling time of muons is still 
shorter than the decay time scale for the relevant energy range. Thus, the conversion efficiency 
into high-energy photons is still large in the Continuous Model, too. 

2.4. Electron-Positron Pair Production 

As was argued in subsection 2.3, a large fraction of the energy of accelerated protons 
is converted into high-energy photons both in the Discrete and Continuous Models. Typical 
energies of photons from neutral pions are above 10 13 eV while those from charged pions and 
muons are below 10 9 eV. The high -energy photons emitted from 7r°s and 7r ± s create electron- 
positron pairs through collisions with residual thermal photons. The cross section of the pair 
creation is written as a± = aig{y) (Berestetskii et al. 1982), where 



(3-y 4 )lni±^-2y(2- 2/ 2 ) 



(15) 



The dimensionless value y is defined by y 2 = 1 — (2m 2 c 4 )/[e 7 e^(l — cos#)]. In figure 6, we plot 
the time scale and the optical depth (r 77 ) against the pair creation in the 100 eV radiation field 
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Fig. 6. Time scale (solid line) and the optical depth (dotted line) of photons for the pair creation in the 
thermal photon field. We assume T = 100 eV and / = 10 9 cm in this figure. 

for the Discrete Model. Figure 6 shows that the photons with energy from 10 8 eV to 

4 x 10 14 eV create electron-positron pairs before they escape from the shell with / = 10 9 cm. 
We define the threshold energy e T , above which photons create pairs. Photons from the decay 
of neutral pions and high-energy part of the synchrotron radiation from charged pions can be 
converted into pairs, while low-energy synchrotron photons will directly escape. 

In the Continuous Model, e T turns out to be similar, since the larger available size 
compensates for the smaller photon density, which results in a similar value of r 77 . Since both 
models give a similar input of high-energy photons, and a similar optical thickness to the pair 
absorption, the subsequent development of the cascade is also similar. 

First we consider pair production of photons from decay of the neutral pions. Such 
electrons and positrons emit the synchrotron radiation the energy of which is given by 



where j± is the Lorentz factor of the pairs. Since the quantum effect is important for electrons 
higher than 10 13 eV for B = 10 6 G, the emitted energy is comparable to the energy of the 
emitter themselves. Those photons create pairs again. Reiterating the pair creation and the 
synchrotron radiation, these high-energy photons cascade into low-energy pairs. 

The energy of electrons/positrons which emit photons with a critical energy of e T = 3 x 10 8 
eV is 8 x 10 10 eV. Thus, the synchrotron photons emitted from electron/positrons above 8 x 10 10 
eV are absorbed through pair production, while those emitted by pairs below 8 x 10 10 eV escape 
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Fig. 7. Photon energy e± in the comoving frame for 772 = 1 and A7 = 1 in Continuous Case. 

from the shell. Thus, for m e c 2 ^± > e T = 3 x 10 8 eV, the injection of pairs due to pair absorption 
of synchrotron photons emitted from higher energy pairs balancing with synchrotron cooling 
makes the resultant number spectrum of pairs as steep as N± oc 7^ 3 . On the other hand, below 
3 x 10 8 eV, pair injection becomes unimportant, and the number spectrum of pairs is determined 
only by radiative cooling as N± oc 7± 2 . Then, the typical energy of the synchrotron emission 
from pairs with break energy e T = 3x 10 8 eV is 



As a result, the photon index of synchrotron radiation becomes —2 above e±, while below 
that the photon index becomes —1.5. Since the energy of photons will be blue-shifted by the 



relativistic motion of the shocked shell [T m ~ r r (T R — y — 1) ~ 200], this typical energy is 
observed as a few hundred keV, which is consistent with the observations. 

In the Continuous Model, similar processes are generated, and the main differences lie 
in the behavior of B and function of r. In figure 7 we plot the typical photon energy 

e± emitted from the pairs of typical energy e T against r. Here, we have assumed r]2 — l, A7 = 1, 
and /b = 0.1. The value e± has been calculated for radii where the photopion production can 
occur by the thermal radiation field. Since e T does not change very much against the change 
in the optical depth, e± practically behaves as oc B oc r" 1 in figure 7. The value e± is around 
0.1-10 keV in the wind frame, which is typically observed as 0.01-1 MeV, similar to the Discrete 
Model. Thus, the results in the Discrete and Continuous Models are basically the same. In both 
models the gamma-ray emission can be ignited by photopion production against the residual 
radiation field. 

Next, we consider the fate of synchrotron photons emitted by charged pions and muons. 




keV. 



(17) 
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Since the energy range is typically 10 5 -10 9 eV, only the high-energy portion of the spectrum is 
absorbed and produce pairs at around e T = 3 x 10 8 eV, while the lower energy photons directly 
escape. Thus, around the break energy, synchrotron radiation of pairs originating from charged 
pions superposes on the smoother component of cascades from neutral pions. 

It is to be noted that this is not the whole story because the thus-produced photons 
overweigh the residual photons and modify the cascading behavior itself. Thus, we expect that 
the optical thickness to the pair absorption drastically increases and that lower energy photons 
which have been regarded to directly escape in the above will also be pair-absorbed to enhance 
the cascade process. 

3. Emission from Cascading Particles 



As discussed in section 2, the residual radiation field of the fireball can trigger a cascade 
from shock-accelerated relativistic protons. The photon field from the cascading particles itself 
will overwhelm the thermal radiation field as the cascade process goes on, and the cascade 
process proceeds nonlinearly. The energy density of the photons generated in the cascades is 
larger than that of the residual photons by about one order of magnitude. In this section we 
discuss the effects of the self-generated photon field on the resultant emission spectra. 

3.1. Self- Generated Photon Field 

Self-generated photons drastically decrease the threshold energy of pair creation (e T ) 
and make inverse Compton scattering as important as the synchrotron radiation as the cooling 
processes of pairs. The cascade process proceeds as shown in figure 8. Although the whole 
process is time dependent, we here treat effects of these photons within the steady state one- 
zone model. In the Discrete Model, for simplicity, we assume that the fraction f 1 of the internal 
energy U- m in the shocked shell goes into the radiation field, and the photon number density 
obeys a power-law distribution as 

fjUml {e^/e y a (for e 7 > e ) 
el | (for £ 7 <e )- 

Of course, this photon spectrum will not generally be coincident with the radiation field that is 
finally realized. As shown below, however, the typical energy of the cascading particles is not 
very sensitive to the shape of the photon spectra, but is determined mainly by / 7 . 

In the Continuous Model we assume that the photon density is similarly obtained using 
U in equation (8). 

3.2. Photopion Production 

For the photon field given by equation (18), r n is analytically obtained as 



n(e 7 ) = (a-2)^\ ^ 7! Z ^ ( 18 ) 
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Fig. 8. Emission process when the photon field due to the cascading particles itself overwhelms the 
residual radiation (see section 3). The notation is the same as in figure 1. 



M7 P ) = G/ 7 2 — TT 2 



a-l 



(19) 



where cr e ff (a) = / dxX~ a<J ir{x) an d 0"cfr(2.2) ~2x 10 _31 cm 2 . The dimensionless factor C for the 
Discrete and Continuous Models becomes T R — 1 and 1Z _1 , respectively. It is noted that the 
expression is independent of the parameters M and A etc., because we have normalized the 
physical values at r = i? p h. 



Since the threshold energy E PiT for which r n = 1 has a dependence of E PtT oc e 



(2-a)/(o-l) 
) 

it turns out to be almost independent of as long as a ~ 2. For a = 2.2 and Eq — 1 keV, all 
protons for which energy is larger than 

-5/6 



E p , T = 10 



13 



0.3 



eV, 



(20) 



produce photopions. The threshold energy becomes an order of magnitude lower than in the 
previous section. Thus, the self-generated power-law photon field enhances the energy transfer 
from accelerated protons to photons and leptons. Since the inelasticity is about 0.2, the energy 
range of produced pions becomes 10 12 -10 16 eV. 

3.3. Photon-Photon Pair Absorption 

The optical depth to the pair creation for a photon with energy e 1 in the photon field 
given by equation (18) is analytically obtained as 
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T ^ ) = Cma) ^(0)°-\ (21) 

where 

F{a) = ^-^ f dyy{l-y 2 Y- 2 g{y) (22) 
a + 1 Jo 

is a dimensionless function [F(2.2) ~2x 10 -2 ]. 

If f 1 ~ 0.3, photons with energy > 1 MeV suffer from pair absorption, and do not directly 
escape from the shell. Thus, the pair cascade fully develops over the whole energy region of 
relativistic pairs. Synchrotron photons emitted from relativistic pairs distribute continuously 
over a wide energy range down to the ~ 0.01 eV without any significant breaks. In the observer 
frame this energy is blue-shifted to a few eV, which is much less than the observed break energy. 
In the Continuous Model, since r 77 ocr -1 , the threshold energy e T increases with r. However, the 
break energy of the synchrotron radiation does not become large enough because the magnetic 
field strength also decreases with r. Even for r = 100_R p h with the fiducial parameters, e T ~ 10 
MeV and e± ~ 0.01 eV in the comoving frame. 

Before we discuss the detailed behavior of the cascade, we examine if photons from pions 
having energies above ~ 1 MeV. 

3.4- Photons from Pions and Muons 

The produced 7r°s promptly decay into two gamma-rays. As is the case of the residual 
thermal photon field, these photons with energy above ~ 10 12 eV cascade into lower energy 
pairs. 

As for the charged pions, high-energy pions cool through Compton and synchrotron 
radiation before decaying into muons, although Compton cooling is fairly suppressed by the 
Klein-Nishina effect in this energy range. In a similar way to that discussed in subsection 2.3, 
we obtain the Lorentz factor of pions, 7^0, whose lifetime is comparable to the cooling time as 
2 = 3m n c C ( mn \ 2 , . 

77K) 4^,0 (/ 7 + / B )m p c 2 UJ ' 1 } 

where the length scale C for the Discrete and Continuous Models becomes Z/(Tr — 1) and 
RpiJ^/Vy respectively We should note that the typical value of £ (~ 10 11 cm) in the Continuous 
Model is much larger than that (~ 10 9 cm) in the Discrete Model. The pions, whose energy 
is larger than 7 7r0 m 7r c 2 ~ 2 x 10 13 (£/10 9 cm) 1//2 [(/ 7 + / B )/0.4]~ 1//2 eV, should cool down before 
they decay into muons and neutrinos. According to equation (14), synchrotron photons from 
charged pions above 3 x 10 13 eV suffer from pair absorption, while those from charged pions 
below 3 x 10 13 eV do not. Thus, a major portion of the energy of high-energy charged pions will 
be injected into the cascade process, and some part results in the emission of photons in the 
sub- MeV range. 

Although charged pions below ~ 10 13 eV decay into muons and neutrinos before sig- 
nificant radiative cooling, the muons take a large fraction (~ m^/m^ ~ 0.8) of the original 
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energy of pions and contribute to photon emission. Muons cool down through the synchrotron 
and inverse Compton mechanisms. As is seen when we estimate 7^0, it turns out that the 
muons, whose energy is larger than "f^m^c 2 ~ 10 12 (£/10 9 cm) 1//2 [(/ 7 + /b)/0.4] _1//2 eV, should 
cool down before they decay into electrons. Since the Klein-Nishina suppression works for 
target photons of energy larger than 10 keV, we should take this limit into account for / 7 . In 
the Discrete Model, since muons are injected in the energy range of ~ 10 12 — 10 13 eV, muons 
lose more than a half of the initial energy radiatively before decay. The cooling time is 
^,cooi — 0.03(£/10 9 cm) 1 / 2 [(/ 7 + /b)/0.4] -1 / 2 s that is comparable to the dynamical time scale 



This energy should be blue-shifted by the relativistic motion of the shocked shell (T m ~ 200) 
and observed as ~ MeV, which is slightly higher than the observed break energy of GRBs. 
The synchrotron spectrum from muon cooling is broadly concentrated around the above typical 
energy. The muons, whose energy is smaller than ^^m^c 2 , decay into electrons. These electrons 
emit very high-energy photons that cannot escape from the shell because of pair absorption. 

In the Continuous Model, since £ ~ 10 11 cm, the typical energy of muon synchrotron 
becomes an order of magnitude larger than in the Discrete Model. 

3.5. Pair Cascade 

In the cascade process, as long as the energy density of power law photons is compara- 
ble to that of the magnetic field, the Compton scattering also produces high-energy photons, 
although Klein-Nishina suppression is significant in the high-energy regime. As the energy of 
pairs in the cascade process decreases, the Compton scattering plays a major role. In addition, 
created pairs increase the optical thickness to the scattering. These processes may modify the 
resultant spectrum. In order to obtain the detailed spectrum in the pair-cascade model we need 
to simulate electron-positron plasmas including the effects of synchrotron radiation, Compton 
scattering, pair creation, pair annihilation, Coulomb scattering, e-e bremsstrahlung, and so on. 
However, such numerical calculation is beyond the scope of this paper. 

Here, we qualitatively discuss the spectrum emitted from cascading pairs, especially the 
possible breaks in the spectrum. We simply estimate the prospective physical conditions of 
the plasmas with the help of published numerical simulations made for pair cascade models for 
compact X-ray sources (Guilbert et al. 1983; Lightman, Zdziarski 1987; Coppi 1992; Stern et 
al. 1995). These previous simulations assumed that copious soft photons are injected at lower 
energies, which correspond to the residual thermal photons of the fireball in the present case. 



l/c. 



The typical energy of muon synchrotron radiation is 




keV. 



(25) 



(24) 
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In the present case, the magnetic field is relatively large; its energy density is comparable to 
that of injected high-energy photons. 

In terms of the compactness parameter, our case corresponds to 
L^a T m p -- f .(Cf±\ 

l °=i^ = Cf ^^ 550 [T3 ' (26) 




where we have taken a volume I because our case is not spherical. Although this typical 
value of the compactness parameter is somewhat large compared with those adopted in the 
above-mentioned studies, the qualitative features will not be different. 

The created pairs above ~ 1 MeV will distribute in a monotonic fashion. The pairs 
cool down mainly through synchrotron radiation and inverse Compton scattering to form a 
power-law spectrum with a number index of around 3, as discussed in subsection 2.4. This 
corresponds to the case when all of the injected energy flows down in the energy space while 
making lower energy pairs and photons. The number spectra of synchrotron emission and 
inverse Compton scattering is described by a power law of index around 2. This very rough 
consideration is a fairly good approximation as the simulations with the power law injection of 
pairs by Coppi (1992) showed the power law indices of photons in the case of / 7 ~ /b become 
around 2 in the range x = e^jvn^c 1 = 10 _3 -10 -1 , which we are interested in. Below this energy, 
the effects of cooled pairs appear while above x — 1 the effects of pair absorption play a role. 

As the low-energy pairs accumulate, a sharp break due to the synchrotron self-absorption 
appears in the spectrum in the case of /b — / 7 . Assuming N± oc 7± < ' P+1 ' ) above E± iVO i a , we obtain 
the synchrotron self-absorption frequency (Rybicki, Lightman 1979) as 



0.4 J \0.3j \0.lj V10 9 cm/ VlMeV/ 
for p = 2.2. The predicted self-absorption energy becomes ~ 2 keV in the observer frame for 
r m ~ 200, lower than the observed break energy. 

Next we discuss the cooled pairs. Following the convention, we define the "pair yield" 
as PY = J d'j±P( r Y±)/ J c?7±Q(7±)(7± — 1), where Q(j±) is the rate of injection of pairs of 7± 
via the pion decay, and P(j±) is the rate at which secondary pairs of j± are created via pair 
production. Since the numerical simulations indicated that for l e ^> 10 the value PY becomes 
~ 0.1 independent of other parameters, we assume PY ~ 0.1 hereafter. This means that 10% 
of the injected energy goes into the rest mass energy of cooled pairs and the remaining 90% 
of energy goes into radiation. Created pairs may cool down to subrelativistic energies and 
thermalize before annihilating. In the steady state, the pair creation rate is assumed to balance 
with the rate of pair annihilation. We thus obtain 
3 If 

-a T cn + n_ = -j P(-?±)dj±, (28) 

where n + and n_ are the number densities of cooled thermal positrons and electrons, respec- 
tively. The Thomson optical depth to thermal pairs is written by 

19 



The equilibrium temperature of the thermal pairs, T±, is mainly determined by requiring that 
no net energy be transferred between the particles and photons. Numerical simulations for 
compact X-ray sources in the past showed G = T±/m c c 2 ~ 10~ 2 -10~ 3 for Z c ^> 10 regardless of 
the other parameters. Therefore, the Compton y parameter, y ~ 4Qt± ^ 1, in our model and 
relatively mild Comptonization may occur. The photon energy tends to be peaked around at 
1-10 keV, which is expected for the observed break energy of GRBs when beaming effects are 
taken into account. 

Because of the large Thomson optical depth of thermal pairs {t±), the mean energy loss 
of photons with energy lower than m e c 2 in one scattering by cold thermal pairs is about e 2 /m e c 2 , 
and the mean number of scattering is about t±. A significant depletion of the spectrum may 
occur in the region 1/t± 1 (Lightman, Zdziarski 1987; Coppi 1992; Stern et al. 1995), due 

to downscattering of photons by cold thermal pairs. The break energy, m c c 2 / t\ ~ a few or a few 
tens of keV, can be the same order of observed break energy owing to the relativistic blue-shift. 
Thus, cooled thermal pairs may possibly produce spectral break compatible with observations 
through thermal Comptonization and downscattering of high-energy photons. However, in order 
to quantitatively investigate this possibility, detailed numerical simulations will be required. 

4. Discussion 

We have shown that Fermi accelerated protons with > 10 13 eV efficiently bring their en- 
ergy into pion production and subsequent photon emission. The photon spectra of the cascades 
from relativistic protons are mainly characterized by synchrotron and inverse Compton com- 
ponents emitted from pairs. In addition, muon synchrotron radiation can make an important 
contribution. Cooled thermal pairs can modify the photon spectra through Comptonization 
and downscattering of high-energy photons. In figure 9, we show the schematic spectra which 
we discussed above. Although we have discussed several mechanisms to reproduce the ob- 
served break energy scale including muon synchrotron radiation, synchrotron self-absorption, 
Comptonization, and down scattering by cooled thermal pairs, at present we have not definitely 
identified any specific one with the observed break feature. The problem remains to be open. 

Let us briefly discuss this problem. The simplest way to produce a sharp spectral break 
is to prevent electrons/positrons from cascading below the typical energy, e T . As discussed 
in section 2, if the photon density in the shocked region is sufficiently low, the pair-creation 
process cannot provide electron-positron pairs under energy e T . The value of e T should be 
large enough (~ 100m c c 2 ) to reproduce the spectral break [see equation (17)]. However, the 
condition that the emission from the cascading particles themselves is more luminous than the 
residual radiation of the fireball implies f 1 ^ 0.1 consistent with the value we have adopted in 
section 3. In addition, / 7 may be larger than ~ 0.1 to reproduce the energetics of GRBs. Thus, 
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Fig. 9. Schematic photon spectrum in the observer frame. The dotted and dashed lines are the com- 
ponents due to the electron-positron pairs and muons, respectively. The solid line shows the integrated 
spectrum. Here, we assume that the energy of the component from the muons is much smaller than the 
energy of the component from the pairs. The break due to the self-absorption appears at ~ 1 keV. The 
downscattering by cold thermal pairs makes both the components above ~ 1 MeV decrease. 

the condition / 7 > 0.1 is universally required so that the pair cascade develops over the whole 
energy region of relativistic pairs (namely e T ~ m e c 2 ). 

As shown in figures 5 and 7, in the Continuous Model the cascade process can occur even 
far outside of i? p h- If internal shocks occur far beyond R p h, the photon density is diluted and 
e T can be large. As discussed in subsection 3.3, however, the smaller magnetic field prevents 
the break energy from increasing sufficiently 

We have neglected effects of the directly accelerated electrons. If the total energy of 
electrons directly accelerated is more than a tenth of the energy of the protons, the synchrotron 
emission from the electrons can make the ignition of photopion production easier. However, 
the emission from the electrons does not modify the cascade process very much. Therefore, 
even in this case e T remains too small to reproduce the observed break energy. 

One may think that the cascade process can be largely altered if /b <C 0.1. Here, we argue 
that this is not the case. As far as B > 100 G with /b >~ 10~ 9 , protons can be accelerated 
beyond 10 13 eV for I = 10 9 cm and photopion production against residual photons occurs 
efficiently. While synchrotron cooling of charged pions is strongly suppressed, neutral pions 
decay into two gamma-rays and create electron-positron pairs through collisions with residual 
thermal photons. The energy of those pairs is comparable to the energy of7r°s (~10 13 eV). Even 
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though the Klein-Nishina suppression makes inverse Compton process ineffective for such pairs, 
the pairs cool down through synchrotron radiation; even for B = 100 G the cooling time of pairs 
via synchrotron radiation is shorter than the dynamical time scale because of their large Lorentz 
factor. The energy of synchrotron photons emitted from the pairs is larger than 100 MeV for 
B > 100 G. Therefore, the synchrotron photons create secondary pairs again (see subsection 
2.4) and the energy of n°s is transfered to lower energy pairs via the cascade process illustrated 
above. The low energy pairs cool through inverse Compton process against residual photons 
and produce the self-generated photon field. When the self-generated photons overwhelm the 
residual photons (/ 7 > 0.01), the cooling processes of low-energy pairs, charged pions, and 
muons are dominated by inverse Compton process, rather than synchrotron radiation. For 
/ 7 > 0.1 the values of 7^0 [see equation (23)], 7^0, and e T remain close to the values for /b = 0.1. 
Thus, even when / B <C 0.1, the cascade process proceeds in a similar way as discussed in section 
3. Cooled muons with 7^0 decay and provide electrons/positrons within a narrow energy range 
around 10 11 eV (see subsection 3.4). Photons emitted from such electrons/positrons can be in a 
narrow energy spectrum, and escape without producing pairs. Even for B = 100 G, the energy 
of the synchrotron photons emitted from the electrons of 10 11 eV is 40 keV in the comoving 
frame. This energy is an order of magnitude larger than the typical break energy. Since the 
Klein-Nishina suppression is significant for such electrons/positrons, the cooling time also tends 
to be longer than the dynamical time scale. In any case, most of the self-generated photons 
come from neutral pions. The spectrum component due to the cascade process initiated from 
7r°-decay may not reproduce the break energy. 

From the above discussion, pairs generally distribute over a wide energy range, which 
overpredicts the flux of low-energy photons between ~ 1 keV and ~ 1 MeV within a simple 
analysis. There may be no physical process to prevent electrons/positrons from cascading into 
lower energy. The possible mechanisms to make the spectral break may be Comptonization 
or downscattering of high-energy photons by cooled thermal pairs. It is difficult to estimate 
the precise effects of Comptonization and downscattering analytically. The Comptonization 
may resolve the problem of the spectral index below the break. We have assumed /b ~ / 7 in 
our discussion. A smaller value of /b makes inverse Compton scattering dominant rather than 
synchrotron radiation. Some numerical simulations showed that the spectral index decreases 
with decreasing /b// 7 (Zdziarski, Lamb 1986; Coppi 1992). A smaller value of /b is favorable 
to decrease the spectral index below the break. After all, we need to simulate the cascade and 
emission processes numerically to predict the final spectra quantitatively. 

Our model may explain some puzzling features of GRBs. For example, some spectra 
obtained by the Ginga satellite show break energies of a few or a few ten keV and spectral 
number index smaller than 2 below the breaks (Strohmayer et al. 1998). These break energies 
are close to the self-absorption frequency in our model. Since the emission spectra in the X-ray 
and gamma-ray range should depend on details of the model parameters, such as the positions of 
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the internal shock and photosphere, our model may provide a wide range of spectral and timing 
behavior. It is interesting to see if our model can reproduce not only the standard spectral 
feature, but also the anomalous spectrum, like the "X-ray rich GRB" (Heise et al. 2001; Kippen 
et al. 2001) or GRB precursor (Koshut et al. 1995). 

Our scenario also predicts several features of neutrino emission. The cooling of muons 
and pions results in a clustering of neutrino energies. From the typical energy of pions at the 
decaying moment, we obtain the energy of which is produced from the decay of ir + 

(tt~) as 0.27 7r0 ^7rC 2 r m ~9x 10 14 (T m /200) eV in the observer frame. If / 7 and /b are sufficiently 
small, the energy of neutrinos can be large. In this case, however, a small f 1 means a weak 
flux of neutrinos. The energy of v e and (p e and z/ M ), which is produced from the decay of /x + 
(/i~), is as small as ~ 7 x 10 13 (T m /200) eV in the observer frame. 

5. Conclusions 

We have investigated the consequences of proton acceleration in the residual thermal 
photon field in the internal shocks of GRBs. Even if the internal energy of heated protons cannot 
be converted into the energy of electrons, radiation energy can be drained from the accelerated 
protons through the photopion production process against residual thermal photons. We have 
examined two extreme models concerning the distribution of residual thermal photons. In the 
Discrete Model, photons are confined within multiple thin shells, while in the Continuous Model, 
photons are emitted from the wind photosphere. In both the Discrete and Continuous Models 
the residual radiation field in the fireball is shown to be able to ignite the photopion production 
and the generated photon field further enhances the photon-initiated cascade. Although the real 
situation lies between these two extreme models, there are no large differences in the ignition 
processes. The Fermi-accelerated protons with > 10 13 eV bring their energy into pions. If the 
number fraction ( of protons are accelerated and distributed as n p oc E~ q , and if they have the 
fraction £ of the internal energy, the energy fraction of pions is obtained by 



The particles cascading from the pions emit photons over a wide energy range. The 
emission properties for the two models are also similar. The spectrum is composed of mainly 
two components: emissions from the electron-positron pairs and muons, produced by high- 
energy photons and charged pions, respectively. The pairs form a power-law spectrum with a 
number index of around 3 over the whole energy region of relativistic pairs. On the other hand 




(30) 



(31) 



for q = 2.2. 
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muons cool down to a characteristic energy before they decay. The radiation due to the pairs is 
composed of a power-law spectrum with an index of about 2 with possible break features around 
the 1-10 keV range in the comoving frame due to Comptonization and downscattering by cooled 
thermal pairs. This may correspond to the observed break when the relativistic beaming effects 
are taken into account. The muons can emit X and gamma-rays by synchrotron process and 
produce a sharp break at around 1-10 MeV in the spectrum. 
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